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We discuss an approximation for the dynamic charge response of nonlinear spin-1/2 Luttinger 
liquids in the limit of small momentum. Besides accounting for the broadening of the charge peak 
due to two-holon excitations, the nonlinearity of the dispersion gives rise to a two-spinon peak, 
which at zero temperature has an asymmetric line shape. At finite temperature the spin peak is 
broadened by diffusion. As an application, we discuss the density and temperature dependence of 
the Coulomb drag resistivity due to long- wavelength scattering between quantum wires. 
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I. INTRODUCTION 

There is theoretical consensus,^ supported by accumu- 
lating experimental evidencepEl that in one dimension 
electrons decay into fractional excitations carrying ei- 
ther spin or charge, called spinon and holon. At low 
energies, interacting one-dimensional (ID) systems are 
described by the Luttinger modelfS which predicts that 
the collective spin and charge modes are decoupled and 
propagate with different velocities. Away from the low- 
energy limit, spin-charge separation holds in the sense 
that spinon and holon branches can still be identified 
in some momentum-resolved experiments such as angle- 
resolved photoemission,Sl but charge and spin degrees of 
freedom are inevitably coupled by dispersion nonlinear- 
ity. A direct consequence is that at finite en ergi es spin 
excitations can contribute to charge responses."™ 

Recently the effects of nonlinear dispersion in dy- 
namical properties of Luttinger liquids have been 
emphasized. S In particular, the interplay of band curva- 
ture and interactions is essential for the interpretation of 
Coulomb drag experiments in parallel quantum wires. 8 9 
As discussed by Pustilnik et al., w the Luttinger model 
cannot account for the leading contribution to the drag 
resistivity when there is a density mismatch between the 
wires, in which case interwire backscattering processes 
are exponentially suppressed at low temperatures. The 
other type of low-energy process, long- wavelength scat- 
tering, is ineffective within the Luttinger model because 
the dynamic charge structure factor (DCSF) S(q, w) for 
small wavevector q is given by a delta function peak at the 
energy of a free boson. In this approximation, the DCSFs 
of two wires with different densities have no overlap and 
the drag resistivity vanishes. For spinless fermions,^ it is 
known that nonlinear dispersion is responsible for broad- 
ening the DCSF i nto a rectangular line shape with width 
proportional to q 2 EES This effect restores a smooth den- 
sity dependence of the drag resistivity. 

In this work we study the DCSF of spin-1/2 fermions 
in the limit of small q and at zero magnetic field. Our 
motivation comes from the search for Luttinger liquid 



behavior in experiments with vertically coupled quan- 
tum wires, in which drag is enhanced by a smaller in- 
terwire separation and the densities in each wire can be 
tuned independently!^ We are interested in the possi- 
bility that spinons give a contribution to the drag resis- 
tivity via spin-charge coupling at finite energies. This 
effect can not be described by the Luttinger model since 
it depends on violating particle-hole symmetry. In order 
to calculate the DCSF, we follow the approach of treat- 
ing band curvature as a perturbation to the Luttinger 
model, and resort to refermionization of the collective 
modes in the cases where perturbation theory is singular. 
While the drag response depends mostly on the spectral 
weight and width of the peaks as a function of wavevec- 
tor and temperature, we also discuss other features of 
the DCSF that are of general interest for the dynam- 
ics of spin-1/2 fermions. These features could be directly 
probed by momentum-resolved techniques, such as Bragg 
spectroscopy in cold Fermi gases £^ We show that at zero 
temperature the charge peak has a q -scaling width, like 
in the spinless case, but there is also a peak due to spin 
excitations which resembles the dynamic spin structure 
factor (DSSF) of Heisenberg spin chains. At zero tem- 
perature, the DCSF diverges at the lower edge of the spin 
peak as a power law with exponent /i s _ = — l/2 + (D(q 2 ). 
At finite temperature, the spin peak is broadened by dif- 
fusion. 



The paper is organized as follows. In Sec. II we present 
the linear response formula for the drag resistivity and 
discuss its relation to the problem of the DCSF of spin- 
1/2 fermions with nonlinear dispersion. In Sec. Ill, we 
derive the effective bosonic Hamiltonian including irrel- 
evant operators associated with band curvature, some of 
which couple charge and spin degrees of freedom. In Sec. 
IV, we describe the line shape of the DCSF in the limit of 
small q at zero temperature. Finite temperature effects 
are also discussed. In Sec. V, the approximation for the 
DCSF is applied to calculate the density and temperature 
dependence of the drag resistivity. Finally, we summarize 
the results in Sec. VI. 
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II. COULOMB DRAG AND DYNAMIC 
CHARGE RESPONSE 

The drag resistivity between two capacitively cou- 
pled wires of length L is defined as the ratio r = 
— {e 2 /2ttK)V2/ I\L, where V 2 is the voltage induced across 
wire 2 (called the drag wire) when a current I\ is driven 
through wire 1 (called the drive wire). For the typical 
setup, see Refs. 18191 Let us assume clean wires (L smaller 
than the mean free path due to impurities) and temper- 
ature regime ksT 3> hvFi/L, where Dpi, i = 1,2, is 
the Fermi velocity for electrons in each wire. The latter 
condition rules out finite size effects which are known to 
produce oscillations in the drag response as a function of 
drive voltageP^ Hereafter we set K= ks = 1- In the lin- 
ear response regime, the drag resistivity at temperature 
T is given bjfffll 



U 2 



dq / dui 
Jo 



q 2 A 1 {q,u))A 2 {q,uj) 
sinh 2 (w/2T) : 



(1) 



where U is the inter wire Coulomb interaction, Vi is the 
charge density and Ai(q,ui) is minus the imaginary part 
of the retarded density-density correlation function in 
wire i — 1,2. Eq. ([I]) expresses the drag resistivity as 
a functional of the dynamic density-density correlation 
function of two decoupled wires. Due to Boltzmann fac- 
tors, the nonzero response comes from the overlap of A\ 
and A 2 integrated up to frequencies of order T. At low 
temperatures compared to the Fermi energies e F i of the 
wires, and neglecting interwire backscattering, the main 
contribution to the integral in Eq. ([TJ) is due to small-g 
(or forward) scatteringP^] 

Our goal is to derive an aproximation for A{q 1 ui) in 
a single wire in the limit q <C kpi and uj <C e F i- From 
this point until Sec. IV, we will be concerned with the 
dynamic response of a single wire and will omit the wire 
index i = 1,2. The wire index will be restored in Sec. V 
when we return to Eq. ([TJ) to compute the drag resistiv- 
ity. 

In order to describe the intrawire interactions, we 
consider a Galilean-invariant model with electron mass 
m and short-range density-density interaction potential 
V(x) 



H 



1 

2m ./o 



dx ¥dl^ 



L pL 

dx / dy V (x — y)n(x)n(y) . 

n JO 



(2) 



Here $ = (ip^,ip±) is a two-component fcrmionic field 
and n(x) — W(x)^(x) is the local charge density. At 
zero magnetic field, the number of electrons with spin 
a =f,l is = N± = N/2. The average density is 
v = N/L and the Fermi wavevector is k F = nv/2. We 
assume the interaction potential to have a finite range R, 
due to screening by nearby gates. For simplicity, we will 
discuss the properties of A(q, uj) in the thermodynamic 
limit. 



The spectral function A(q,uj) in Eq. ([TJ) satisfies the 
fluctuation-dissipation theorem 



2A(q,u) = (1 



,-u/T 



)S(q,uj), 



(3) 



where S(q, ui) is the DCSF given by 



S(q,u)= dxe- iqx dt e lut (n(x,t)n(0,0)) . (4) 

Jo J -oc 

Since S(—q,u) = S(q,u), hereafter we take q > 0. Since 
we are interested in the regime T <§C cf, we shall start by 
discussing the DCSF at T = 0. The exact result for the 
nonintcracting case, V — 0, is 



S (q,oj) = (2m /q) 6(q 2 /2m - \u) - v F q\), 



(5) 



where vf — kpjra is the Fermi velocity. For V = 0, 
the spectral weight vanishes outside the particle-hole 
continuum defined by the upper and lower thresholds 
UJ ±(q) — v fQ i q 2 /2m. As a function of energy u, the 
line shape of the DCSF in this case consists of one rect- 
angular peak whose width is given by the band curvature 
scale 5u)(q) — q 2 /m. 

It is important to note that for spin- 1/2 fermions the 
DCSF at small q cannot be obtained by perturbation the- 
ory in the interaction. In fact, the first order correction 
to S(q,uj) has logarithmic divergences atdi« &±(q) 



SS(q,u) m(2V q -V ) 



S (q,v F q) 



In 



nq 



UJ + 



(6) 



where Vk is the Fourier transform of V(x). These di- 
vergences signal edge singularities and are reminiscent 
of the result for spinless fermionsJiSl but there is an im- 
portant difference. For spinless fermions, the prefac- 
tor of the logarithmic divergence is m(Vo — V q )/irq; for 
short-range interactions, Vq — V q ~ q 2 , it vanishes as 
q — > 0. In contrast, the prefactor in the spinful case, 
m(2V q — V~o)/nq s» mVo/nq, diverges as q — > 0. The dif- 
ference stems from the amplitude for s-wave scattering 
between electrons with opposite spin. This indicates that 
the limits q —¥ and Vq in S(q, uj) do not commute. 
In the limit q -C mVo, it is important to treat interac- 
tions exactly and account for the effects of spin-charge 
separation, as we will discuss in the following sections. 



III. EFFECTIVE MODEL FOR SPIN-CHARGE 
COUPLING 

In the regime q -C mVo, we can treat interactions ex- 
actly and regard band curvature as a pert urbation , with 
q/k F playing the role of a small parameter! 12 * 17 * 1 ^ In the 
linear dispersion approximation, bosonizatiorP of Hamil- 
tonian ^ is standard and leads to the Luttinger model 
with Hamiltonian density 



2ttv c (J 2 + J 2 ) 



2ttv s 



2Trv s gJ R ■ 3 L 
(7) 
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Here J r /l (Jr/l) are chiral U(l) charge (SU(2) spin) 
currents, v c (v s ) is the charge (spin) velocity, and g is 
the bare coupling constant of the marginally irrelevant 
backscattering operator. The long- wavelength part of the 
charge density fluctuation is 

n = 2^K~ c (J R + J L ), (8) 

where K c is the Luttinger parameter for the charge sec- 
tor. Galilean invariance implies K c — vp/v c . At weak 
coupling, Vq -C Vf, we have v c » Vp + Vq/it, v s s» vf 
and g ~ V-nep/irVF-io- semiconductor quantum wires, 
typical values of K c ss 0.7 have been reported™ 

The charge and spin currents can be expressed in terms 
of chiral bosonic fields as 



Jr/l = Td x <p c R/L /V^, 
Jr/l = Td x <p R/L /V^, 



(9) 
(10) 



which obey the commutation relations 

^fiS(x — x'). The transverse part of the spin currents 
can be written as 



J 



R/L 



= J 



1 



R/L 



j- _ jx 
J R/L — J R/L 



IJr / l ~ 2ir 6 

■TV - JL 

R / L ~ 2tt 6 



+iy/ilt<p s R/L 



(11) 

(12) 



where the short-distance cutoff is set to 1. 

The leading (dimension-three) perturbations to model 
([7]) , generated by the quadratic term in the electron dis- 
persion as well as irrelevant interactions, are 



SH = (4^ 2 /3) [r,_(J% + Jl) - V+ (J 2 rJl + J 2 l Jr) 
+k-(J r 3 2 r + J L Jl) + k+(J r J 2 l + J L J 2 R ) 
+K3(Jl + Jr)3l-Jr]. (13) 



The last three terms in Eq. ( |13[ ) couple spin and charge. 
Importantly, only even powers of the spin currents are 
allowed in 8% due to SU(2) symmetry. Direct bosoniza- 
tion of Hamiltonian ^ produces all terms in Eq. (13 1, 
except for the K3 term. This does not mean that K3 
vanishes (such term is allowed by symmetry), but rather 
that it must be generated at second order in the electron- 
electron interaction. 

In fact, we can derive phenomenological relations for 
all coupling constants. The exact parameters rj± can be 
related to the change of v c and K c under a shift of chemi- 
cal potential \i. The calculation is analogous to the spin- 
less case in Ref. H21 simplifying for the case of Galilean 
invariance where K c — vp/v c and using the result for the 
compressibility dv/dji = 2K c /irv c , we find 



17- 



2VK, 
3 (\_ 
m 



2JK r 



1 dv c 

1" V c~^~ 

m 0[i 



(14) 
(15) 



Likewise, an infinitesimal chemical potential shift Sfi 
modifies the spin velocity v s by giving a finite expec- 
tation value to the charge currents (Jj) = (J R ) = 
6n\/K c /(2irv c ) in the k± terms in Eq. (13 1: 



47r 2 2ir 

— («_ + k+)(J r )(3% + J|) = T ^ S (J| + J 2 L ). (16) 



This relation fixes the sum 

K- + k+ = (vj \f~K~c) dvjdix. 



(17) 



Moreover, Galilean invariance requires that the charge 
current and momentum operators for model ^ plus ( 13 1 
be proportional to each otherPThe momentum operator 
is obtained from the energy-momentum tensor; its den- 
sity is 



V 



2k F 

~7k, 



:{Jr-Jl) + MJ 2 r-JI) 



2n 



(J|-J|). (18) 



The current density SJ{x) is obtained from the continuity 
equation for the charge density: 

d t n{x) =-i I dx' [n(x),Hi(x')+SH(x')] = -d x J(x). 



(19) 

We then impose the condition J(x) = V(x)/m for a 
Galilean-invariant system. The relation between the co- 
efficients of the spin contributions to V{x) and J{x) leads 
to 

(20) 



Eqs. (17 1 and (20 1 allow one to determine n± by simply 



measuring the spin dispersion at low energies. 

Finally, the coefficient K3 is related to the variation of 
the backscattering coupling constant v s g under a change 
of the chemical potential, 



K3 



3v c d(v s g) 
2^TC C dn 



(21) 



Since g is marginally irrelevant, we expect K3 to be more 
irrelevant than the other coupling constants in Eq. (13 I, 



in the sense of logarithmic corrections to scaling. This 
will be discussed in the following subsection. 



Renormalization group flow with irrelevant 
operators 



All the operators in Eq. (13 1 are irrelevant and have 
the same scaling dimension x — 3. The renormalization 
group (RG) equations for the irrelevant coupling con- 
stants (including the marginal g in Eq. ^) can be de- 
rived by integrating out high-energy modes in the par- 
tition function as one lowers the ultraviolet momentum 
cutoff A. Following the notation of Ref. 19, we define 
"dimcnsionless" coupling constants (which have dimen- 
sions of velocity) fj± = Ar/±, R± — Ak± and £3 = AK3. 
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To obtain the quantum corrections to scaling, we use the 
operator product expansion (OPE) of the spin currents 1 



:J a L (z)::J b L (0) 
: J£(z)::4(0) 



6 ab 


i£ abc 


8n 2 z 2 


2nz 


S ab 


i£ abc 


8ir 2 z 2 


h 2-kz 



Jr(0) (22) 



where z — v s T+ix and z — v s r—ix, with r the imaginary 
time, and e abc is the Levi-Civita antisymmetric tensor. 
[The normal ordering symbol :: is implicit in the Hamil- 
tonian Eqs. Q and The OPE for the charge fields 

is simply 



: J l (x,t) :: J L (0,0) 
: J r (x,t)::J r (P,0) 



1 



8tt 2 (v c t + ix) 2 
1 

87T 2 (v c T — ix) 2 



(23) 



We integrate out high-energy modes in the shell 1/ A < 
\z\ < 1/A' with A' = Ae~ M , dl < 1. This choice of 
cutoff is rotationally invariant for the spin modes, but 
elliptical for the charge modes. In order to get a nonzero 
contribution to the RG equation after integrating out the 
shell in the (x,t) plane, it is important to contract both 
right and left movers for a given species (charge or spin) 
at the same time. 

In the presence of the dimension-three operators, the 
velocities, Luttinger parameter and chemical potential 
are renormalized, but flow to their fixed-point values in 
the low-energy limit. This flow is already taken into ac- 
count if we use the exact parameters. The interesting RG 
flow here is given by the coupled equations for k 3 and g 



dg 
dl 
dk 3 

~dJ 



-(l + 2g)k 3 . 



(24) 
(25) 



There are no corrections to the scaling of fj± and k± to 
second order in the coupling constants. Eq. (251 can be 
rewritten as 



dn 3 
di 



-2ffK 3 . 



(26) 



On the right-hand side of Eqs. (24) and (26 1 we have 



terms of zeroth order in A. It follows that 

d In g 1 d In k 3 

dl ~ ~ 9 = 2 dl ' 

The solution is of the form 

^((■)/W)\ 2 = const.. 



(27) 



(28) 



The scaling of the marginal coupling constant is the fa- 
miliar one 



5(A) 



l+.gln(A /A)- 



(29) 



where Ao is the initial value of the cutoff. As a result, 
for positive g <C 1 and for A <C A^er 1 / 9 the effective 
g(A) vanishes logarithmically as g(A) ~ l/ln(A /A). If 
the bare g at high energies is of order 1, the perturbative 
result in Eq. (29 1 is still valid if g is interpreted as g(Ao) 
at some scale Ao <§C kp such that g(Ao) -C 1. In any case, 
we obtain g(A) ~ l/ln(A /A) in the low-energy limit. 
More interestingly, Eq. (p8| implies 



«s(A) = 



«3 



[l + 9 ln(A /A)P 



(30) 



Therefore, k 3 (A) vanishes as n 3 (A) ~ l/ln 2 (A /A) as 
A — > 0. This will be important in Section |IVD| when 
we compare leading logarithmic corrections to the DCSF 
due to g and n 3 . 



IV. DYNAMICAL CHARGE STRUCTURE 
FACTOR 



by 



In the bosonized form of Eq. 



the DCSF is given 



S{q,u) = -8K c lm C rct (q,tu) 



(31) 



where C lct (q,ui) is the retarded correlation function for 
the charge current Jr + Jl, which can be obtained by an- 
alytic continuation from the Matsubara correlation func- 
tion 



C(q,iw) 



C a/3 (q,iuj) 



C a p(q,iu>) 



E 

,/3=R/L 
fL 

dx e- iqx 

'0 JO 

x{J a (x,r)Jp(0,0)). 







dre 1 



(32) 



(33) 



Eq. (33) involves the charge boson propagator. Within 
the Luttinger model, the charge boson is free and we have 
C° LL = C L , C° RR = C R and C° LR = C° RL = 0, with 

C r / l (x,t) = (Jr/l(x,t)J r / l (0,0)) 

07t z [v c t =p ixy 
Taking the Fourier transform, we obtain 

1 ±q 



47T iuj =p v c q 



(35) 



As a result, the DCSF calculated in the linear dispersion 
approximation is given by 



S(q,u) = 2K c qd(uj - v c q). 



(36) 



That the DCSF is given by a delta-function peak at 
the energy of the free charge boson follows from spin- 
charge separation and Lorentz invariance of the Luttinger 
model. This should be contrasted with the free-electron 
result in Eq. §5§ , where the peak associated with particle- 
hole excitations has a q 2 broadening due to the curvature 
of the dispersion about the Fermi points. 
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A. Width of the charge peak 

We can calculate S(q, ui) beyond the Luttinger liquid 
result by analyzing the effects of the boson-boson inter- 
actions in Eq. (13). First, let us consider the broaden- 
ing of the charge peak. The charge-only rj± terms are 
familiar from the spinless caseP^ They account for the 
decay of one charge boson into two charge bosons. In 
particular, r/_ is a three-leg vertex in which one right- 
(left-) moving boson decays into two right-(left-)moving 
bosons, thus coupling the single-boson state to degen- 
erate multi-boson states. It is known that perturba- 
tion theory in fy_ is ba dly divergent, but can be dealt 
with by refermionizationl^^ Near the charge mass shell, 
uj k, v c q, we introduce a spinless holon field iprR such 
that R tp c R = \f2J R . The ?7_ term in Eq. jl3j) then 
maps onto a parabolic dispersion about the holon Fermi 
point 



47T 2 

IT 



V-JR 



v- 
2V2 



€r9 2 A 



c,R- 



(37) 



It can be argued that 7/_ determines the exact broad- 
ening of the DCSF to order q 2 because it gives rise to 
the bosonic diagrams that are most singular at u) = 
v c q^ Within the approximation of neglecting the other 
dimension-three operators, the charge sector of the Lut- 
tinger model plus the r\_ term refermionizes into a free 
fermion model with dispersion e c (k) ss v c k+rj_k 2 / (2^/2), 
for k measured from the right Fermi point. The support 
of the charge peak in the DCSF is then given by the spec- 
trum of excitations with a single holon-anti-holon pair. 
Due to the curvature of the dispersion, these excitations 
define a continuum bounded by 



u c ±{q) = v c q ± T7_g 2 /2v / 2- 



(38) 



Therefore, at order q 2 , the charge peak h free- 
fermion-like line shape 



2V2K C Q { V _q 2 

S{q,w) ~ -—7= - \w- v c q\ 

V-Q V 2V2 



(39) 



The parameter y/2/r)_ can be interpreted as a renormal- 
ized holon mass. 

It is interesting that the limits q — > and Vq —> in 
the width do not commute. For q <C mVo <C kF, we have 



from Eq. (14 1 that 5w c (q) = r^q 2 /V^ « q 2 /(V2m). The 



\/2 factor makes the charge peak narrower than the free 
electron result in Eq. (|5])j231ln particular, this means that 
in the regime q <C mVo the holon dispersion (which shows 
up, for instance, in the single-electron spectral function) 
should not be regarded as a smooth continuation of the 
electron dispersion. An important crossover happens at 
q ~ mV . 



The result in Eq. ( 38 I can be directly compared with 



the exact width of the two-holon continuum for an in- 
tegrable model, such as the Yang-Gaudin modelPH We 




ja 

•J n 



ja 

■J n 



Figure 1: Decay of a charge boson (propagators denoted by 
wiggly lines) into two right-moving spin bosons (propagators 
denoted by straight lines). This process leads to a spin peak 
in the dynamic charge structure factor. 



have numerically solved the standard Bethe ansatz inte- 
gral equations for the spectrum of elementary excitations 
of the Yang-Gaudin model. We verified that the width 
defined as the difference between the maximum and min- 
imum exact energies of a particle-hole excitation in the 
holon Fer mi sea for momentum q <C m(v c — v s ) agrees 
with Eq. ( 38 I , including the factor of \[2 and with cal- 



culated from the phcnomcnological relation in Eq. ( 14 1 



Corrections to Eq. ( 39 I due to residual holon- holon in- 



teractions are higher order in q. These corrections include 
a high fre quency tail at order 77^, analogous to the spin- 
less caseji^^UJy anc l possible asymmetries of the charge 
peak near the edges of the two-holon continuum, due to 
x-ray edge type singularities.^ But before we discuss the 
behavior near cu c ±(q), we turn to the contributions from 
the spin operators in Eq. (13). 



B. Spin peak 



The k± operators in Eq. (13 1 allow for decay of the 



charge boson into two spin bosons moving in the same 
direction and carrying the total energy uj sa v s q. The 
corresponding three-leg vertices are illustrated in Fig. [T| 
As noted in Ref. [THJ this process gives rise to a narrow 
peak in the DCSF centered at u> = v s q which corresponds 
to a charge-carrying spin singlet excitation. Let us cal- 
culate the correction to the charge boson propagator in 
Eq. (33 1 using second order perturbation theory in k±. 
For instance, the 0(fti) correction to Crr is 

8G K RR {q, iw) = 32^ V [C R (q, iuj)} 2 IL RR (q, iu), (40) 

where we made use of the identity : J 2 R / L ■= 3(Jj^/ L ) 2 
and introduced the boson self-energy 



U a/3 (q,iuj) = - J dxe lqx 



x / dre 1 ^ S a (x,T)S p (x,T). (41) 
Jo 

Here S a , a — R, L, are the free chiral spin boson propa- 
gators 

S ah S R/L (x,r) = (J a R/L (x,T)J b R/L (0,0)) 
S ab 1 

8TT 2 (v s t =p ix) 2 ' ^ 
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In momentum and frequency space, we have 

1 ±q 



Sa./L{<l,iu) 



An iu) =F v s q 



(43) 



We then calculate Hrr that appears in Eq. ( 40 ) 

dq' [°° dJ q' 



(47r) 2 n fliJ ( ?! ia;) = - 



2tt 



2ir iu' — v s q' 



icu — ioj' — v s q + v s q' 
7 3 1 



127T iui — v s q 



(44) 



Adding up all second order contributions from k_ and k + 
and taking the imaginary part of the retarded self-energy, 
we obtain^ 

8S(q,u) « [K c /12){a- + a+) 2 q 3 S(u J - v s q), (45) 
where 



Ot± = K±/(v c ± V s ). 



(46) 



This shows that the DCSF exhibits a narrow peak with 
spectral weight of order q 3 at the spin mass shell u> = v s q. 
The question then is how the spin peak in Eq. ( 45 1 is 



broadened by treating band curvature operators to higher 
orders. We first note that the condition q <C m(v c — v s ) 
ensures that the spin peak is well separated from the 
charge peak. By analogy with the discussion in Sec. 
|IV A[ we expect that the width of the spin peak is set 
by decay processes that couple degenerate states with 
multiple spin bosons propagating in the same direction. 
However, in contrast with the case of the charge peak, 
in Eq. (13 1 there is no dimension-three spin-only chiral 



operator that would be equivalent to a parabolic disper- 
sion about spinon Fermi points. Thus the broadening of 
the spin peak must come from higher-order on-shell de- 
cay processes. In fact, the leading irrelevant spin-boson 
interactions allowed by symmetry are quartic in the spin 
currents 



C-(4/l) 2 ,Ch 



R ' Jl) 2 ; A 2 J_R ' Ji(J^ 



•Ji) 



(47) 

In principle, these operators are present as perturbations 
to model plus Q. They are also generated by "pro- 
jecting" into a subspace with energy \lu— v s q\ -C (v c — v s )q 
and integrating out "high-energy" charge bosons. The 
spin part of the resulting model for uj w v s q is equivalent 
to the low-energy effective model for the XXZ spin chain 
at zero magnetic fieldPS 

Pe rtur bation theory in the dimension-four operators in 
Eq. (47) is highly singular for ui ~ w s 9- Unfortunately, 



it is not known how to sum up the expansion in this case. 
Refermionization does not solve the problem because the 
effective fermionic model with dimension-four operators 
contains not only band curvature terms, like ^ d 3 i]j , but 
also intrabranch (i.e. which do not mix R and L) resid- 
ual interactions of the form i/j^dxip dxip^ip, which also 



contribute to the broadening at leading order in q. Nev- 
ertheless, simple power counting tells us that the width 
of the spin peak should scale like 5uj s (q) ~ 0(q 3 ), rather 
than 0(q 2 ). This is consistent with the result for the 
DSSF of the XXZ model at zero field where it is known 
that the spectral weight is dominated by two-spinon ex- 
citations and the exact spinon dispersion takes the form 
e s (k) = v a sin(fc) v s (k — fc 3 /6 + . . . ) about the spinon 
Fermi points. 



C. Edge singularities of the spin peak 

To be able to say more about the line shape of the 
spin peak, we refermionize the spin currents into inter- 
acting spinless fermions. This is equivalent to inverting 
the Jordan- Wigner transformation in the continuum and 
writing down a SU(2) symmetric model for the fermions 
associated with the spin excitations. In other words, the 
idea is analogous to deriving the bosonic Hamiltonian 
for the Heisenberg spin chain by starting from the XXZ 
model and tuning the Luttingcr parameter to the SU(2) 
symmetric value (a strongly interacting limit with Lut- 
tinger parameter K — 1/2), as opposed to deriving the 
bosonic Hamiltonian directly from the Hubbard model 
(in which case the spin bosons come out noninteracting 
with K s = 1) .111221 The new ingredient here is that the 
spin excitations are coupled to gapless charge modes. 
The mapping of the bare chiral fermion densities to 



the spin currents in Eq. (j^FJ) is ^ 



R/L^s ,R/L 



(3J 



R/L 



J£, iJ )/2, as follows from a canonical transformation for 
the spin bosonic fields. Spin inversion symmetry implies 
that the dispersion of these fermionic spinons is particle- 
hole symmetric. We assume that the exact dispersion 
about the right Fermi point is given by 



,(*) 



(48) 



with the unknown parameter 7 > 0. We expect that 
7 stems from dimension-four operators in the bosonic 
model and is of order \/(rnkp). 

In terms of fermions, the operator 3 2 R that gives rise to 
the spin peak in Eq. (45) creates particle-hole pairs on 
the spinon Fermi sea. We can study the behavior near 
the edges of multi-spinon continua using the methods of 
Refs. I16|24l The absolute lower threshold o; s _(g) = e s (q) 
is defined by an excitation with a particle at the Fermi 
surface and a hole at momentum —q below the Fermi 
point. For u> — uj s -(q) jq 3 , we define a "deep spinon" 
subband by expanding^ 



(49) 



where %j}\ r creates low-energy spinons near the right Fermi 
point and d\ creates a hole at momentum — q below the 
Fermi point. This leads to the quantum impurity model 



ri a 



U e + d\(u s . 
+2Trq(K,'_J r 



- iud x )d s - (V r J* 
K,' + Ji)d\d s , 



ViJ?)d\d s 
(50) 
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where u ~ v s — 3jq 2 is the velocity of the d s hole. The 
spin-only part of the quantum impurity model given by 
the first line in Eq. ( 50 I is derived as explained in Ref. |2"H 



by applying the mode expansion ( 49 1 to a generic model 



of interacting spinless fermions with the dispersion in Eq. 
( 48 1 . The second line stems from the coupling of the en- 



ergy density of the spinon field to the bosonized holon 
density. Here J r /i stand for the bosonized charge cur- 
rents with cutoff at energy scale <C jq 3 , which allows the 
"high-energy" spinon to emit low-energy charge bosons 
such that the energy remains near oj — ui s -(q). Note 
that in this procedure we keep only marginal operators 
in the quantum impurity model, as irrelevant operators 
can only introduce subleading power-law singularities at 
the threshold. This is not to be confused with the pres- 
ence of irrelevant operators in the original bosonic model 
( 13 1 , which are essential to argue for the nonlinearity of 



holon and spinon dispersions and for the very existence 
of the deep spinon threshold. 

Rather than keep track of the parameters in the deriva- 
of model ( 50 ) , it is more useful to introduce the 



tion 



model phenomenologically (it contains all the marginal 
operators allowed by symmetry) and to fix the coupling 
constants by symmetry and phenomenological relations. 
The parameters V r n can be fixed by realizing that the 



same model (50 1 can be used to calculate the lower 
edge singularity of the DSSF. This is because the op- 
erator JJw^ that enters the longitudinal spin-spin corre- 
lation function also creates two-spinon excitations in the 
fermionic representation, and the lower threshold of the 
support of the DSSF is also given by the deep spinon 
excitation. The important constraint comes from SU(2) 
symmetry, which imposes that the exponents for the lon- 
gitudinal and transverse DSSF must coincide!^ We re- 
produce this argument in detail in the appendix. 

The coupling constants k' ± in Eq. ( 50 I are related 



to exact phase shifts at the holon Fermi points due to 
the creation of a high-energy spinon. We want to show 
that at small q these are also related to the band curva- 
ture parameters in Eq. (13). It is easy to show that an 



infinitesimal change in the chemical potential 5[i gives 
rise to a shift in the energy of the high energy spinon 
6lo s - = 5^l\/K c (k'_ + n' + )q/v c . This allows us to write 



duj. 



fW c q dpi 



(51) 



But from the exact spinon dispersion we have lo s 
v s q + 0(q 3 ), hence 



dv. 



+ 0(q 2 ). 



(52) 



The second relation for be obtained by im- 

posing Galilean invariance to Hamiltonian ( 50 1 . Simi- 
larly to the discussion in Section III, we compare mo- 
mentum and current operators. The contribution from 
the spin-charge coupling terms in Eq. ( 50 1 to the cur- 



rent density (defined from the continuity equation for 



the charge density) is 



Jd = \fK c q{fi!_ -K' + )d\d s . 



(53) 



Therefore, if we consider an excited state in which we 
create a particle-hole pair of spinons with a deep hole 
at momentum kp — q and a particle at kp, the current 
of this state is y/K c q(n'_ — k' + ). Demanding that this 
current be equal to the momentum q of the state divided 
by the mass m, we find 



k' + = l/( v / i^m). 



(54) 



Comparing Eqs. ([52j and (|54| with (|17() and (|20|, we 

conclude that 



«4 = k± +G(q 2 ). 



(55) 



Using model ( 50 I and Eq. ( 55 I , we can show (see ap- 



pendix) that the DCSF diverges at the lower edge of the 
two-spinon continuum as S(q,uj) ~ (lu — w s _) A ' s -, with 
exponent 



fi s _ = -1/2 + (a 2 _ + a 2 + )q 2 /2 + 0(q 4 ), 



(56) 



where a± is defined in Eq. (46). Therefore, as q — > 0, the 



exponent approaches the universal value —1/2, which de- 
pends only on SU(2) symmetry. The q 2 correction to n s - 
is due to the coupling to gapless charge bosons with en- 
ergy <C 7<7 3 . This exponent should be contrasted with the 
square-root singularity of the DSSF for the Heisenberg 
modelPO We note that /i s _ differs from the correspond- 
ing exponent for SU(2) bosons at the magnon threshold. 

Mm = _l + 0( ? 2 ) .M 

The upper edge of the two-spinon continuum is given 
by uj s+ (q) = 2e s (q/2). As discussed in Ref. [24j near 
this edge the spectral weight is suppressed by resonant 
scattering between spinons with equal velocity. If most of 
the spectral weight of the spin peak is due to two-spinon 
excitations, the width can be defined as 



5ui s (q) = io s+ (q) - u> s -(q) = 37g 3 /4. 



(57) 



While the upper threshold of two-spinon continuum in 
the integrable XXZ model exhibits a square-root cusp, 
here we expect that the upper threshold of the two-spinon 
continuum is rounded by higher order (in q) processes, at 
least for non-integrable models. 



D. Smearing of the charge peak: decay rate of the 
charge boson 

In general, we expect S(q, w) to have nonzero spectral 
weight anywhere above the lower threshold td s -(q). A 
tail between the spin and charge peaks is generated due 
to the decay of the charge boson into a pair of L and R 
spin bosons, as depicted in Fig. [2] The effective vertex is 
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Figure 2: Processes that contribute to the effective vertex k| 
in Eq. ( |59[ ) to leading order in dimension-three operators. 
The decay of the charge boson into a pair of right and left- 
moving spin bosons leads to a tail in the charge structure 
factor for uj > v s q. 



calculated from the three-point function 

G(k l5 k 2 ,k 3 ) = ((J L + J i? )(k 1 )J£(k 2 )4(k 3 )) 
4tt 2 

= ^-Kf8 ab (C L + C fl )(k 1 )^ i (k 2 ) 
x^(k 3 )(2^) 2 ,5(k 1 +k 2 + k3), (5f 



where k = (k,uj) is a two-momentum. To first order 
in K±,ft3 (leading order in q/kp in the contribution to 
S(q,(jj)), there are two contributions to the effective ver- 
tex, one from K3 and the other from a combination of k± 
and g. We find 



«3 



-M K - 



(59) 



Away from spin and charge mass shells, i.e. for \u> — 
v s q\ ^ Suj s (q) and \us — v c q\ ^> 5ui c {q), the tail of the 
spin peak can be calculated by second order pertur bation 
theory in Kg ff . Similarly to the calculation in Sec. 



IV B 



we obtain a correction to the charge boson propagators 

16tt 4 



SC K 4 (q, 



{K?yC a {q,iu) 
xCp(q,iu;)U RL (q,iu}), (60) 



where a, ft — R, L and Hrl (q, iu>) is the self-energy with 
one right-mo ving and one left-moving spin boson, as de- 
fined in Eq. (41 ). The calculation of Hrl yields^ 



n&(g,w) 



1 



32tt 3 



2v* 



2 2 

v* q A 



8w? 



x lot 



(v s q) 2 - (ui + irff 
4v?A 2 



(61) 



While the real part is ultraviolet divergent , the imaginary 
part is not. The imaginary part gives the tail in the 
DCSF 



24u? 



c ff\2 r 



v c q 



-1 2 



v 2 q 2 



(co 2 ~V 2 q 2 ), (62) 



for ui > v s q and \ui — v c q\ ^> r/_q 2 . 

At this point, we recall that both g and K3 that ap- 
pear in the amplitude for scale logarithmically with 



the infrared cutoff (see Sec. Ill A I. Within RG improved 



perturbation theory, the bare g and K3 in Eq. (59) are 
replaced by the renormalized ones in Eqs. ( pfj| and (30 1, 
with cutoff set by the small momentum A ~ q. In the 
long-wavelength limit, such that g(q) ~ 1/ \n(kp /q) <C 1, 
we have 



g{q) ~ l/Hk F /q) 
«3(?) ~ (k 3 / g 2 )/ In 2 (k F /q). 



(63) 
(64) 



At leading logarithmic order, w e ca n drop the contribu- 



tion from k 3 in ^(q) and Eq. (62) becomes 



6S{q,u) 



3K c [g(q)} 2 
32w? 



v 2 q 2 



I .2 2 2\ 

(uj -v s q ) 



(65) 

The DCSF was calculated by similar methods in Ref. 
[TBI but the tail between the spin and charge peaks was 
not obtained because backscattering processes (g and K3 
in our notation) were neglected. 

The presence of the tail means that the charge peak 
discussed in Sec. |IV A| is inside a continuum of spin exci- 
tations. The coupling to the continuum results in a decay 
rate for the charge excitations. The imaginary part of the 
self-energy Hrl can be absorbed into the charge boson 
propagator in the form 



±q 



4-7T uj =p v c q + i/t c 



where r c 1 (q) is the decay rate given by 



1 



n(nf) 2 (v 2 - v 2 )q 3 
192w? 



In the limit g(q) ~ 1/ h\.(kp/q) -C 1, we obtain 

1 = 3[g(q)} 2 n(^+ K+ ) 2 (v 2 -v 2 )q 3 
t c 256^3 



(66) 



(67) 



(68) 



For q -C m(v c — v s ), the decay rate in Eq. (68 I is smaller 
than Su> c (q) ~ q 2 , which is due to decay of the charge 
boson within the charge sector (see Sec. IV A I. However, 
r^ 1 is important because it is responsible for rounding off 
the edges of the charge peak. This can be confirmed by 
calculating the decay rate for a single high-energy holon - 
a d c particle in the quantum impurity model for the edges 
of the two-holon continuum, similar to the calculation in 
Ref. [30] for the spinless case.. The decay rate is due to 
the perturbation 



Sn 3 ~ £_ 3 dld c d x if R d x 



<Pl- 



(69) 



The parameter £ 3 gives the amplitude for a process in 
which a holon scatters off two spinons moving in opposite 
direction. This is not a three-electron scattering process 
and in principle £ 3 ^ even for integrable models. Using 
the result in Eq. (5.8) of Ref. 1301 we obtain 



1 



oc 



(& 2 (v 2 c-V 2 ){A(qW 



(70) 



9 



where A(q) is the cutoff of the high-energy subband. Set- 
ting A(q) ~ q. we recover the momentum and interaction 
dependence of the decay rate in Eq. ( 68 1 if we assume 



that £3 does not vanish as a power law of q in the limit 
q — > 0. From comparison with Eq. (67 1, we expect 



3 OC K3 



if 



This should be contrasted with the spinless 
case, where the coupling constant in Eq. (5.13) of Rcf. 
[3"U] has to vanish like q 2 because of statistics, since there 
is no s-wave scattering for spinless fermions. For the 
spinful case, statistics alone does not imply that the am- 
plitude £3 vanishes as q 2 or high powers of q. 



Note also that the result for 1 /r c in Eq. ( 68 1 is non- 



perturbative in the electron-electron interaction, since for 
q <C mVo <C kF, we have r" 1 ~ g 2 Voq 3 /kp, which is third 
order in the interaction strength. Furthermore, this re- 
sult implies that, even for an integrable model, the power- 
law singularities 16 at ui c ± are removed at order q 3 . The 
same decay rate l/r c rounds off the singularity at the 
holon mass shell in the electron spectral function This 
is remarkably different from the spinless case, where it is 
believed that integrable models can have exact singular- 
ities above the lower threshold because the decay rate of 
a high-energy particle may vanish exactly.^! 

We note that if the phcnomcnological relations 



Eqs. (17 1 and (21 1 are valid for the running coupling 
constants and we substitute them in Eq. ( 59 1 , we 



find nf = -{3/2)v c v s K~ 1/2 dg/d^ At low energies, 
5(A) ~ g(A )/[l+g(A ) ln(A /A)]. It is not clear how im- 
posing the phenomenological relations at all energy scales 
can be reconciled with the result from the RG. It is im- 
portant for our results in Eqs. (65) and (68 1 that even 
if dg/d/i — at some energy scale, the effective vertex 
k| will be generated by the RG flow because K3 and g 
scale differently, and the leading logarithmic dependence 
is due to the g term in Eq. ( 59 1 . 



Finally, putting together all the pieces, we construct 
the final picture for the DCSF at zero temperature in 
Fig. [3] Note, in particular, that there is only a rounded 
threshold at u) ~ lu, — 



E. Finite temperature effects 

We now discuss the effects of finite temperature on the 
broadening of the charge and spin peaks in the DCSF. 
This will be important to compute the temperature de- 
pendence of the drag resistivity in Eq. (|T|) . 

We consider the regime where both temperature and 
band curvature energy scale are small compared to the 
scale of spin charge separation: T <C (v c — v s )v c /ri- and 
q <C (v c — v s )/r)—. Assuming that v c , v s . and v c — v s are 
all of the order of vf and 1/rj— is of the order of m, these 
conditions mean roughly T <C Tp = mv 2 F /2 and q^kp- 
This is the regime in which we may expect the spin and 
charge peaks to remain well separated. Neglecting the 
overlap between the spin and charge peaks, the line shape 
of the charge peak can be approximated by the finite 
temperature result for the imaginary part of the density- 



S{q,u>) 



I 
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s+ 



LU 



Figure 3: (Color online) Schematic excitation spectrum (inset 
on the upper right corner) and line shape (main figure) for the 
DCSF of spin-1/2 fermions at small q. The dashed lines in the 
(q, uj) plane represent the linear dispersion of charge and spin 
modes in the Luttinger model. The two-holon and two-spinon 
continua are represented by blue and yellow regions bounded 
by u) c ±(q) and u) a ±(q), respectively. At zero temperature and 
in the regime q <C m(v c — v s ), the line shape of S(q, w) (solid 
red line in the main figure) has well defined peaks inside the 
two-holon and two-spinon continua. At finite temperature 
5u} s (q) <C g 2 T <C v c q, diffusion broadens the spin peak into a 
Lorentzian (dashed green line). 



density correlation function for free fermions with mass 



V2K C 
V-Q 



[n F {w + ) - n F {w-)}, (71) 



where rip (w) = 1/(1 + e u / T ) is the Fermi-Dirac distribu- 
tion function and w± = [(lu ± Slu c /2) — (v c q) 2 ]/(2Soj c ), 
with 5lo c = Sui c (q) = r\-q 2 j\[2. The width of the 
charge peak at finite temperature is then of the order 
of max{77_(7 2 , r]-qT/v c }. 

The calculation of the width of the spin peak is more 
complicated because we do not have an approximation in 
terms of noninteracting spinless fermions. For the pur- 
pose of calculating the drag resistivity in Eq. , we are 
only interested in whether for fixed small q and at low 
temperature the spin peak can become broader than the 
charge peak. On the one hand, the broadening due solely 
to band curvature must be of order ('jq 2 /v s )T. As long 
as q,T/v s <C m(v c — v s ), '7-/7, this is small compared 
to the broadening of charge peak. On the other hand, 
thermal effects have a stronger effect on sp in exc itations 
because the latter are damped by diffusion! 32 ! 33 ! 

Recall that the spin peak stems from the self-energy 
with two spin bosons propagating in the same direction 
(see Fig. 0. We can calculate the finite temperature 
broadening by neglecting band curvature operators and 
applying perturbation theory in the marginally irrelevant 
operator g in Eq. ^ . As we did in section IV D 



we ne- 



glect the K3 vertex in the leading logarithmic approxima- 
tion. To order (n±g) 2 , there are two types of diagrams 
in the self-energy for the charge boson, as illustrated in 
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Fig. [4] The first type amounts to a self-energy correc- 
tion to the spin boson propagator. The transverse part 
of the perturbation, — irv s g( J~£ J~j^ + h.c.) gives rise to a 
nonzero imaginary part of the retarded self-energy, which 
can be calculated following Ref. [M] We can sum up the 
series for this type of diagram by defining the dressed 
spin propagator 



S R (q,iw) = 



1 



(72) 



4tt iuj — v s q — £(<?, iuj, T) 

The other type of diagram (Fig. [3Jb)) is a vertex cor- 
rection. Since the thermal broadening of the spin peak 
is already obtained within the approximation of keeping 
only self-energy-type diagrams like the one in Fig. [I|a), 
we will make the approximation of neglecting vertex cor- 
rections. By doing this, the two-spin-boson correlation 
function becomes 



{4Tr) 2 n RR (q,iuj) 



~ I < ^r T ^S R {q',w n ) 

xSr(q- q' ,iu - iv n ) 

dq' q'(q - q') 



2ir iuj - 



x[n B (v s q') 



- v s q - lu q , 
n B {v s q' - v s q)],(73) 

where n B (u) = l/(e w / T - 1) and E 9 = E,(T) = 
S(q r , v s q,T). The decay rate for the spin boson at fi- 
nite temperature is the well-known spin current relax- 
ation rate^l 



— ^ = -Im E™* 

Ts(T) 



MT)] 2 T, 



(74) 



where g(T) « g/[l + g\n(T F /T)}, with T F ~ mv 2 F , is 
the coupling constant at scale T. The finite temperature 
result for the spin peak in this approximation is then 



SA{q,u,T) 



Kr 



q 3 F(q,T) Ts 



48tt 1 + [(u - v s q)T s /2f 



(75) 



where 



/oo 
duu{\ - u) 
-OO 

x[n B {v s qu) -n B (v s q(u- 1))], (76) 

such that F(q, T — > 0) = 1. Since we neglected band cur- 
vature effects in the spin boson propagator, this approx- 
imation is only valid for 1/t s (T) 'yq 3 . In Eq. (73 I, 
we also assumed l/r s (T) <C v s q. Therefore we expect 
that the line shape of the spin peak at finite temperature 
be well described by a Lorentzian for a range of q that 
scales linearly with temperature, q ~ T/v s , such that 
T/g(T) « K 3 / 7 ) 1/2 . 

More generally, for fixed q there is a crossover temper- 
ature T* given by the condition 1 /t s (T* ) ~ jq 3 at which 
the line shape of the spin peak goes from highly asym- 
metric (with a peak near the zero temperature threshold) 
below T* to approximately Lorentzian in the diffusion- 
dominated regime above T*. 




(b) 




Figure 4: Diagrams at 0(K±g 2 ) in the self-energy for the 
charge boson: (a) self-energy correction to the spin boson 
propagator; (b) vertex correction. The bubble with mul- 
tiple lines denotes the correlation function for the operator 
i^Cv'k-'el), following Ref. [34] 



V. APPLICATION TO THE DRAG 
RESISTIVITY 

Finally, as an application of our results for the DCSF, 
in this section we will discuss the density and tempera- 
ture dependence of the drag resistivity in Eq. 0. Here 
we restore the index i = 1, 2 for properties of the drive 
wire and drag wire, respectively. Note that Eq. (JlJ is 
symmetric under the exchange of drive and drag wires. 
We may then assume v c \ < v C 2- We envision an exper- 
iment in which the electron density of the drag wire is 
varied while the one in the drive wire is kept fixed. In 
principle, this could be achieved with vertically coupled 
wires whose densities are tuned independently by top and 
back gates. We can think that the sharply peaked DCSF 
of the drive wire is the reference for integrating Eq. ([lj 
in the (q, lj) plane, and approximate 



U 2 



\'K 3 V\ViT 



/ dq 2 - 

Jo sinh [v c i 



(v cl q/2T) 



< / dwA 1 (q,u)A 2 (q,w). 
'0 



(77) 



Clearly, the forward scattering contribution to the drag 
resistivity is maximum for ideal density matching, v c \ = 
v C 2^ m which case it is dominated by the overlap of 
the charge peaks of the two wires over modes with lo 
v ci q T <C Tp. The temperature dependence in this 



case is the same as for spinless fermions 
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T 2 for 



m \ v ci ~ v c2 \ < ^ < Tf . Also like the result for spinless 
fermions, a much weaker response, r ~ T 5 , is obtained for 
general density mismatch in the regime T <g. rn\v 2 cl — v 2 2 \, 
from the overlap of the charge peak for one wire with the 
tails of the DCSF for the other wire. 

The new effect due to the spin degree of freedom is 
related to the presence of a spin peak in the DCSF il- 
lustrated in Fig. [3] It suggests that the drag resistivity 
can be enhanced when the electron densities in the wires 
are rather different but the charge peak of one wire over- 
laps with the spin peak of the other wire. This happens 
over an extended region in the (q, oS) plane if v c \ ps d s2 - 
The interpretation is that under this condition holons in 
wire 1 can efficiently scatter off spinons in wire 2, which 
in their turn transfer momentum to holons in the same 
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Let us discuss the temperature dependence of the drag 
resistivity when v c i — v s2 . We are interested in the 
regime T -C mv c \(v c2 — v c i), where spin charge sepa- 
ration may be measurable. According to the result in 
section IV E the spin peak SA(q,uj,T) for values of q 
that scale linearly with T assume a Lorentzian line shape 
at low temperatures such that T/g 2 (T) <C {v^/lz) 1 ^ 2 ■ 
The width of the charge peak for wire 1 for q ~ T/v S 2 
is of order r/^^T 2 /v s2 v c \. If the temperature is also 
low enough that T/[g 2 (T)] 2 <C u c i^s2/?7-,l ~ T F , the 
spin peak for wire 2 is broader than the charge peak 
for wire 1. Using Eq. (71 1 for wire 1 and Eq. (75) 
for wire 2, we fi nd that in this regime the drag resistiv- 
ity in Eq. fn\ scales like r ~ T 5 /[g 2 (T)} 2 . However, 



if the value of g 2 {T) is small even at intermediate tem- 
peratures, there will be in general a temperature regime 
92(Vg 2 /j 2 ) 1 / 2 «T« mv c i(v c2 — v c i) where diffusion is 
not effective, in the sense that for q ~ T/v s2 the broad- 
ening of the spin peak due to diffusion is smaller than the 
one due to band curvature. In this case, the spin peak 
for wire 2 is narrower than the charge peak for wire 1 
for the same value of q <C k Fi . Therefore, the integral in 
Eq. ( 77 1 can be evaluated by considering that the entire 



spectral weight of the spin peak is inside the charge peak. 
In this regime, the drag resistivity in Eq. (77 1 scales like 
r ~ T 4 . In summary, 



0.05 




0.9 0.95 

Figure 5: Drag resistivity as a function of the spin ve- 
locity v S 2 in the drag wire for three values of temperature: 
T/(k 2 F1 /m) = 0.03,0.05,0.07 (bottom to top). In this graph 
we set Vo/(nvFi) = 0.4 and V2k F / (nv F i) = 0.1 in the 



weak coupling expressions for v c 



V- 



for 



wires i = 1,2, using equations in section III. The value of 
r(v a 2) is normalized by the value at zero density mismatch, 
ro = r(v S 2 = v s i), for each temperature. Notice the peak in 
r when the spin velocity of the drag wire matches the charge 
velocity of the drive wire, v S 2 « v i- 
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Notice that in the low temperature limit diffusion sup- 
presses the drag by making the spin peak broader than 
the charge peak. Nevertheless, even in the spin diffusion 
regime the smallness of g 2 (T) ~ l/ln(Tp/T) makes the 
drag larger than the background contribution from the 
tails of the DCSF, which scales like r ~ T 5 , as discussed 
above. 

In order to compute the expression in Eq. ( 77 1 , we 



have used the same strength of the electron-electron in- 
teraction Vq and V 2 k F for both wires. We estimated 
the parameters of the DCSF in each wire using the phe- 
nomenological relations in section III, expanding to first 
order in the interaction. The finite temperature Ai(q,co) 
are approximated by the sum of a dominant charge peak 
given by Eq. (71 1 and a smaller spin peak given by the 
Lorentzian in Eq. (75 1. For modes with q ~ Tjv c \ that 
contribute to the integral in Eq. ( 77 1 but are not in 



the regime jq 3 <C 1 /t s (T) , the Lorentzian is not a good 
approximation for the line shape of the spin peak (which 
must become asymmetric with a peak near the lower edge 
of the two spinon spectrum) . However, the drag resistiv- 
ity is not sensitive to the detailed line shape since in this 
regime the spin peak is much narrower than the charge 
peak for the same q and what matters most is whether 
the spin peak for the drag wire falls inside the charge 
peak for the drive wire. 

Fig. [5] illustrates the dependence of the drag resis- 



tivity on the density mismatch between the two wires, 
parametrized by the ratio v s2 /v c i. The effect of spinon- 
assisted Coulomb drag is observed as a small peak in 
the drag resistivity when the wire densities are such that 
v s2 ~ v c \. The height of the peak relative to the dom- 
inant response at zero density mismatch (y s2 ~ v s \) in- 
creases with increasing temperature. This is because, as 
temperature increases, modes with larger q, for which 
the spin peak in the DCSF have relatively larger spectral 
weight, start to contribute to the integral in Eq. ( 77 1 . On 



the other hand, increasing temperature also broadens the 
dominant peak observed at v s2 ~ v s i, which eventually 
obscures the smaller peak at v s2 v c \. Therefore, it 
seems that the contribution of spinons to Coulomb drag 
would most likely be observed as a shoulder in the den- 
sity dependence of the drag resistivity at intermediate 
temperatures. 



VI. CONCLUSIONS 

We have studied the dynamic charge response for spin- 
1/2 fermions in one dimension beyond the usual linear 
dispersion approximation of Luttinger liquid theory. Un- 
like the spinless case, the limits of small momentum q and 
weak electron-electron interactions Vq do not commute. 
This is due to the interplay of spin charge separation 
and band curvature effects. The problem of calculating 
the dynamical charge structure factor for electrons with 
mass m in the regime q vtiVq cannot be approached by 
perturbation theory in the electron-electron interactions. 
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We have used a bosonized Hamiltonian and discussed the 
effects of irrelevant perturbations associated with band 
curvature. We determined phenomenological relations 
for the coupling constants of these perturbations, includ- 
ing the ones that couple charge and spin dynamics. The 
renormalization group equations for the irrelevant oper- 
ator denoted by k 3 , which couples charge and spin and 
mixes right and left-moving spin modes, shows that its 
effective coupling constant has a nontrivial logarithmic 
scaling in the low energy limit. 
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Appendix A: Calculation of the lower edge exponent 



Based on a picture in which collective charge and spin 
modes can be refcrmionized into spinless fermions (holons 
and spinons, respectively) with nonlinear dispersion, we 
presented an approximate line shape for the dynamic 
charge structure factor valid in the long wavelength limit . 
As a function of frequency, the dynamic charge struc- 
ture factor has a dominant charge peak associated with 
two-holon excitations, whose width scales like q 2 . How- 
ever, the spectral weight extends down to a lower thresh- 
old described as a two-spinon excitation. We calculated 
the exponent of the power-law singularity at this lower 
threshold and found that in the limit q — > it converges 
to the universal value fi s - = — 1/2 + 0(q 2 ), which de- 
pends only on spin SU(2) symmetry. We expect that the 
spectral weight near this lower threshold is largest within 
a two-spinon continuum, giving rise to a spin peak in the 
charge structure factor. There is also a tail of the spin 
peak above the two-spinon continuum. The coupling be- 
tween spin and charge in the vicinity of the charge mass 
shell gives rise to a decay rate for the holon at order q 3 . 
This decay rate is responsible for rounding off singular- 
ities at intermediate thresholds such as the edges of the 
two-holon continuum, regardless of the integrability of 
the model. 

At finite temperature, an important difference between 
charge and spin excitations is that the latter are damped 
by diffusion. In the dynamic charge structure factor this 
effect is manifest in the Lorentzian broadening of the spin 
peak in the regime where the spin current decay rate is 
large compared to the band curvature energy scale for 
the spinons. 

These results allowed us to calculate the Coulomb drag 
response between two quantum wires, taking into account 
the spin degree of freedom. In comparison with the re- 
sult for spinless fermions studied in Ref. 1101 there is an 
additional effect (spinon-assisted Coulomb drag) due to 
spin-charge coupling: at low temperatures the drag resis- 
tivity as a function of density mismatch has a peak when 
the charge velocity of one wire matches the spin velocity 
of the other. The temperature dependence of this drag 
peak has signatures of spin diffusion. 

Note: After this work had been submitted, Ref. [36] 
appeared with related results for the exponent of edge 
singularities at arbitrary momenta. 



We provide details for the calculation of the expo- 
nent of the power law singularity in the DCSF at u> = 
ui s -(q) at zero temperature. While the asymptotic value 
fi s -(q — > 0) = —1/2 is universal and depends only on 
SU(2) symmetry, following the argument of Ref. [SSJ we 
also obtain the q 2 correction due to spin-charge coupling. 

As mentioned in Sec. |IV C| the lower edge of the spec- 
trum is the same for the DCSF and for the DSSFs, defined 



r-h r+oo 

S ab (q,U)) = J dxe~ lqx 1 " - 



dte iut (S a (x,t)S b (0,0)). 

(Al) 

Here S(x) = ^ (x)^^(x), with r the vector of Pauli ma- 
trices, is the spin density operator. Its long wavelength 
components are represented by S = Jl + J^. For fixed 
q, the lower threshold below which the DSSFs (both lon- 
gitudinal and transverse) vanish is controlled by a deep 
spinon excitation with energy uj s -(q) = e s (g). 

The components of S satisfy SU(2) commutation rela- 
tion 



[S*(x), S b {x')} = ie abc S c (x)S(x - x 1 ) 



(A2) 



We will refcrmionizc the spin density to spinless fermions. 
This can be done using an inverse Jordan- Wigner trans- 
formation in the continuum 



S*{x) = ^l(x)^ s (x) + const. 
S + (x) = 9\(x)e 



%-k dx 0(x — x')S z (a/) 



(A3) 
(A4) 



where ^ s is a spinless fermionic spinon field, which obeys 
anticommutation relations {^f s (x),^l(x')} = 5(x — x'), 
and 0(x) is the left-continuous Heaviside step function 
with 6»(0) = 0. 

A generic spin Hamiltonian that is a function of the 
local spin density and respects SU(2) symmetry takes 
the form 



H, = CiS • S + C 2 d x S ■ d x S 



(A5) 



where . . . stands for higher-order irrelevant operators. 
By means of Eqs. ( A3 1 and ( A4 1 , this maps onto a model 



of interacting spinons. The precise form of the spinon 
Hamiltonian is not essential here, but it must be such 
that in the low energy limit it yields the same equal-time 
correlation functions as the spin part of the Luttinger 
model Q. This is directly accomplished if model ^ 



13 



is recovered by bosonization of *f? s , analogously to the 
bosonization of the XXZ modelP^ In this approach, the 
SU(2) symmetric model corresponds to strong interac- 
tions between spinons. Here we assume that the effective 
model for the spinons in the metallic case can be ap- 
proached in the same way as the effective model for the 
Heisenberg spin chain, namely by starting from a general- 
ized model of weakly interacting spinless fcrmions where 
we can expand the dispersion about the Fermi points to 
bosonize the low-energy degrees of freedom. SU(2) sym- 
metry is only imposed at the end, on the results for the 
spin-spin correlation functions, to fix the parameters of 
the effective model. (It is conceivable that the strength 
of the spinon-spinon interaction could be tuned in a mi- 
croscopic model for a metal with spin U(l) symmetry 
and that such weakly interacting limit could be realized.) 
In addition to the interactions in the effective spinon 
model, we must account for the coupling to gapless charge 
modes. For the purpose of deriving the exponent at the 
spinon edge, the latter can be described by bosonic fields 
Pr/l a ^ an steps and need not be refermionized. The 
spin-charge coupling is then equivalent to spin-phonon 
coupling in spin chains. 

Let us consider that in the ground state the spinons 
form a Fermi sea with a particle-hole symmetric band 
(due to spin inversion symmetry S z — ¥ —S z ). The ele- 
mentary S z = excitations are particle-hole pairs in the 
spinon Fermi sea. Similarly, there are particle-hole ex- 
citations in the holon Fermi sea, but in our low-energy 
effective model for the lower edge these are treated as 
charge bosons (since we can neglect band curvature for 
the holons). This picture is supported by the Bethe 
ansatz solution of the Yang-Gaudin or Hubbard models. 
We assume that we can start from a model of noninteract- 
ing spinons with dispersion f< f R / L {k) ~ ±(va k— r y^k 3 + 
. . . ) about the Fermi points ±fcp. That the spinon Fermi 
wave vector is given by kp s — kp follows from the peri- 
odicity of the spin excitation spectrum, which is gapless 
at momentum 2fcj?pSI Another interpretation is that in 
this approach the number of spinons, 



dx ¥ s (x)V s (x) = k Fs L/TT = N/2, (A6) 



is the spinon Hamiltonian with spinon-spinon interac- 
tions contained in H" lt , H c is the charge Hamiltonian 
given by the charge part of Eq. and T-L cs con- 

tains spinon-holon interactions. Due to spinon particle- 
hole symmetry, the latter can only contain irrelevant 
operators (dimension three and higher), for instance, 

^l,R/L d ^sM/L d ^R/L- 

The parameters vf ^ and 7^ are renormalized by in- 
teractions (both spinon-spinon and spinon-holon). We 
denote the parameters of the exact spinon dispersion by 
v s and 7. As usually done for v s , the exact 7 can be 
extracted from the Bethe ansatz solution in the case of 
integrable models. For repulsive electron-electron inter- 
actions, we expect v s < v c , where v c is the exact charge 
velocity. In this case, the lower edge of the spectrum 
of any dynamical correlation function at small momen- 
tum q (such that |7g 2 | <C v s < v c ) is controlled by the 
spinon branch line, with a single "deep spinon" with en- 
ergy e s (q) = e s fl(<z) ~ v s q — jq 3 and a certain number of 
spinon or holon excitations at the Fermi points. 

Whether or not the band curvature of the spinon dis- 
persion can be neglected in the calculation of dynamical 
quantities depends on the frequency range of interest.^ 
Far enough from the threshold, for \u> — e s (o)\ 3> 7<Z 3 , 
we are allowed to drop the band curvature operators 
and the dynamics is captured by Luttinger liquid the- 
ory. We can bosonize ip s ,R/L m the standard wayj^l 
using Tp S: R/L ~ g-^v^^fj/L J ^/2na, with chiral bosonic 
fields 4> S R / L and a a short distance cutoff. The quadratic 
Hamiltonian in terms of spin bosons is diagonalized by a 
transformation to (p s R , L 



± <p' R = 



± , 



'r)i 



(A9) 



where K is the Luttinger parameter for the spinons. The 
latter is fixed by SU(2) symmetry. In Abelian bosoniza- 
tion, we write S z ~ d x (p R — d x (p s L , which has scaling 
dimension 1. For S + , we bosonize the fermion and string 
operators in Eq. ( A4 1 and obtain 



5+ 



)<Pr. 



(A10) 



is fixed by the condition that the state constructed by 
adding (removing) N/2 spinons to the ground state is a 
fully polarized state with N spins up (N spins down), 
which no more spinons can be added to (removed from). 

Expanding the spinon field about ±kp, fy s ~ 
e tkFX rj) sR + e~~ lkFX ip S L, we can write a phenomenologi- 
cal Hamiltonian density of the form 



rl — H s + H c + T~Lc 



(A7) 



where 



n s = 



+1 pt L (zv^d x - ii {0) dl)4> sL + K nt (A8) 



which has scaling dimension K + 1/AK. Demanding that 
this dimension is equal to 1 gives K = 1/2. 

On the other hand, for \uj — e s (q)\ <C jq 3 , the behavior 
of dynamical correlation functions is sensitive to the band 
curvature energy scale. In this regime, it is important 
to consider that spinons created near the threshold with 
energy e s (q) travel with a different velocity than spinons 
at the spinon Fermi surface. It has become standard 
to treat this type of problem using quantum impurity 
models in analogy with the x-ray edge singularity!^ In 
order to describe the deep spinon threshold, we expand 
the spinon field with low and high energy subbands: 



ij sR ~ ^ sr + e-^cf sl + e iqx d t2 . 



(All) 
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The subband momentum cutoff is taken to be small com- 
pared to q. For q > 0, we consider deep holes with mo- 
mentum about kp — q and high-energy particles with mo- 
mentum about kp + q. These two types of excitations are 
degenerate as a consequence of particle-hole symmetry. 
The low energy -0 sr and ip s i fields are then bosonized, 
while d s i t 2 are treated as mobile impurities. We denote 
by ip s r H the chiral spin boson fields with the reduced cutoff 
at scale q. The interactions among low and high energy 
spinons and low energy charge bosons can be described 
by the effective Hamiltonian density 



n = n c + h.j + n d + n sd + n cd . 



(A12) 



Here 1-L C is the free charge boson Hamiltonian with re- 
duced cutoff at scale q (fields denoted by y'r/z) - I n the 
spin-only part of H, 



^ = f[(«) 2 + (a^?) 2 ] 



(A13) 



is the Luttinger model for low energy spin excitations 
(here written in Abelian bosonization form) , 

U d = 4iM«) - iud x ]d sl + (1 -> 2) (A14) 

is the kinetic energy of the high energy spinons. with 
u = u(q) the corresponding exact velocity, and 



n sd = -^(Vtdstf - v r d xip s r )(dl 2 d s 

V47T 



tf Bl d Bl ) (A15) 



is the coupling between high energy spinons and low en- 
ergy spin bosons. The amplitudes V r and Vi scale ~ q 2 
because they represent short-range inter action s between 
spinless fermions and vanish as q — > 0) 10 l 24 l Note that 
H s d is invariant under the particle-hole transformation 

^s,R/L -> ^J,fl/Z> Wllich takeS K/l "> -V'/f d sl d s2- 

Here we have neglected the backscattering operator g in 
Eq. (|7J, which becomes marginal at the SU(2) point; 
however, we point out that it might be important for 
logarithmic corrections to edge singularities, which are 
known to exist for the DSSF of the Heisenberg spin 
chain The spin-charge coupling in H is given by 



with 



V, 



ll/r 



l/r 



V, 



2(v s ±uY ll/r v c ±U r,.±r. 
This transformation takes 



l/r 



V, 



l/r 



(A18) 



d s i- 



d s i,2e 



M =F- 



(A19) 



where d s i.2 are free (up to irrelevant operators). The 7's 
are interpreted as phase shifts at the spinon and holon 
Fermi points due to the creation of a high-energy spinon. 

The exponents for the edge singularities are then calcu- 
lated using the methods of Refs. 16 24] The threshold for 
the longitudinal DSSF S zz (q,oj) is given by the correla- 
tion function for the operator that creates a particle-hole 
pair of spinons with a hole at kp and a particle at kp + q 
(or equivalently a hole at kp — q and a particle at kp) 



with 



d s 2r^sr 



e -tv^f(A r ^*+Aivf+ArV?+ A ?V) e ) ) (A20) 



A, 



and 



V/z 











V2 












V2 


n) 




1 


lr/l 



2yjK c 7T 



(A21) 



(A22) 



Using Eq. (A2CM, we calculate the Fourier transform 
of the correlation function (B z (x,t)B\(Q,Q)) and find a 
power-law singularity S zz {q,uj) ~ (u> — with ex- 

ponent 

Mzz = -l + (A r ) 2 + (A0 2 + (A^) 2 + (An 2 

2 t /1 \ 2 



4+ 1 3 -> 

2 \2 7T 



1 
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2" 



(A23) 



n cd 



v r c d x <p c r )(dl l( i al + dl 2 d a2 ). 



The amplitudes Vfn in Eq. (|A16|) stem from opera- 



tors like ^siidxipsRdxtPr/V taking the high-energy mode 
in the expansion of ipR. As a result, scale like 



q. The relation to the parameters in Eq. (50 1 is 



V, 



l/r - ±y / 2K c irqK , ± . 
Hamiltonian % can be diagonalized by a unitary trans- 
formation of the form H = UHU^ where U = U1U2 with 



The last term amounts to an orthogonality catastrophe 
contribution to the exponent due to coupling of the d s i,2 
particles to gapless charge modes. 

Now consider the transverse DSSF S A (q, to). In this 
case the operator Eq. ( A4 1 creates a particle and a string 



B\{x) ~ 4 e -V-/2«-«) 

~ J^g-iv^fA^.+Ajipf+A^+A^) ^24) 



with 
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The Fourier transform of the correlation function 
(B+(x,t)Bl(0,0)} leads to S+-(q,u) ~ (lj - lu s -Y+- 
with the exponent 



-1 
1 



2 \ 2 
n 















V2 " 





(A26) 



SU(2) symmetry implies that \i zz = ji^ , but this is only 

one equation. However, we can actually get infinitely 
many equations by imposing SU(2) symmetry for the 
singularities that differ from the above by zero-energy 
excitations in which spinons are transferred between the 
Fermi points. These are created by electron backscatter- 
ing processes in the language of spin- 1/2 fermions, which 
are umklapp processes for spinons, 



— i y/An n ( ip f — (p ® ) 



(A27) 



where n is an integer (n < on the right-hand side of 
Eq. (A27) corresponds to the Hermitian conjugate of the 



left-hand side) . The excitations that differ in momentum 
by 2nkp have thresholds at the same frequency u> s -(q) 
because the spin spectrum is periodic in momentum with 
period 2k p. The exponents for \q— 2nkp\ <C kp are given 

by 



f-^zz,n 
— ,n 



-l + (A r ,„) 2 + (A,„) 2 + (A^) 2 + (An 2 , 
-1 + (A' ) 2 + (AL) 2 + (A^) 2 + (An 2 , 



where 



A, v 

K, 



i 

71 
i 

71 
i 

71 
i 

71 



-2n 



-2n 



-2n 



-2n 



(A28) 



The condition \i xz n = (J>-\ , n is satisfied for all n if 

and only if 7 r /7r = ji/ir = 1/2. With this result, the 
exponent in the lower edge of the DSSFs (or the DCSF) 
for q <C kp becomes 



1 



With 7; c /r given in Eq. (Al8l and Vf. = ±V2K~ c irqK'± 



(A29) 



fixed as explained in section |1V C| we obtain the final 
result in Eq. ( 56 I . 
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